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conformation may significantly affect the conductance of these devices. Future 
improvements in device reproducibility and yield will require control of these factors. 
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1. Introduction 

Electron transfer is of fundamental importance to proteins whose functionality relies on 
redox reactions with other proteins, cofactors, or the external environment. In proteins 
containing heme groups, a single Fe ion in a protoporphyrin ring changes valence state 
between ferrous (Fe 2+ ) and ferric (Fe 3+ ) oxidation states as an electron is accepted by 
or donated to the protein [JJ. Much experimental work has been performed to develop 
an understanding of protein electron transfer mechanisms (2J E] most of it based on 
statistical averages of protein ensembles measured by spectroscopic and electrochemical 
techniques (3J E] . 

Electron transfer through individual molecules has been probed recently using 
electrochemical scanning tunneling microscopy (EC-STM) [HE]. In EC-STM, a gate 
voltage is used to align a molecules redox level to the Fermi levels of the tip and 
substrate with the sample in solution at room temperature. The results are explained 
using resonant tunneling [H |6] and two-step electron transfer models [5] . According to 
resonant tunneling models, electrons transfer via resonant tunneling due to alignment 
of molecular orbital energy levels with the Fermi level of the contact electrodes. In the 
two-step process, conformational changes in the molecule occur during tunneling which 
change the orbital energy levels. This process is thought to be especially prevalent 
in redox molecules [HI E]- In a related group of experiments, measurements on large 
assemblies of localized proteins using mercury contacts have been used to distinguish 
between tunneling and hopping electron transfer mechanisms [8]. 

Nanometer-gap electrodes fabricated by electromigration techniques [9] have made 
it possible to study electron conductance through small single redox molecules. By 
applying a bias voltage Vb between electrodes and a gate voltage Vq to the sample's 
substrate, molecular energy levels can be probed and characteristic single-electron 
transistor (SET) behavior, such as a Coulomb blockade or a Kondo resonance, can 
be observed [TQl, [Til [El EES]. m particular, these experiments have demonstrated 
that molecular junctions fabricated in this way can be used to probe the electron 
transfer properties of the molecules. In these measurements, the presence of Coulomb 
blockade phenomena by itself is not proof that the molecular properties are being probed 
because the characteristics of the junction itself after the breaking process can affect the 
measurements [H] , and therefore control experiments must be performed to validate the 
data. These measurements, if performed correctly, rely on a resonant tunneling process 
which in principle indicate that tunneling is a possible electron transfer mechanism. 
Our aim in this work is to determine whether tunneling is a viable electron transfer 
mechanism in redox proteins such as myoglobin by attempting to fabricate and measure 
single electron transistors. 

Here we describe the study of single-electron transistors composed of the protein 
myoglobin (Mb). Mb is a roughly spherical globular metalloprotein containing a heme 
group surrounded by a folded single-chain polypeptide consisting of 153 amino acids 
(apomyoglobin) [15] . Mb was chosen because 1) it is relatively small (4.4 nm x4.4 nm x 
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2.5 nm), which is desirable to increase electron tunneling probability, 2) the apo form 
is readily available and well studied and can thus be used as a control to discern effects 
related to the heme group, and 3) its structure and vibrational spectra are well known. 
In other words, Mb is a model system with which to test bioelectronic devices based on 
heme-containing proteins. 

This paper is organized as follows. After briefly discussing the general properties 
of myoglobin and single-electron transistors in the Background section, we describe 
the experimental techniques used, including break junction fabrication, electrical 
measurements, and protein incorporation. In the Results section we describe results 
from bare junctions, junctions with apomyoglobin (apoMb), and junctions with Mb. In 
the Discussion section we analyze our data using the known properties of myoglobin 
and the results from the bare and apoMb samples. Because SET conductance can 
result from metal grains produced by the electromigration process [T6l IT3] . an important 
distinction between samples with and without protein incorporation is the observation 
of inelastic tunneling from discrete vibrational modes [HI HH [17]. We also identify 
heme-group resonant tunneling by comparing the behaviors of apoMb samples with Mb 
samples. The Conclusions section summarizes our results and identifies areas where the 
experiments can be improved. In Appendix A we describe differential SET conductance 
calculations with and without tunneling-induced structural changes in the protein. 

2. Experimental techniques 

2.1. Break- junction fabrication 

The junctions were grown on Si substrates with a SiC>2 layer 200 nm thick. The resistivity 
of the degenerately-doped silicon wafers ranged from 0.002 fi-cm to 0.003 fi-cm. The 
Si substrate served as the gate contact. The widths of junctions were usually between 
100 nm and 300 nm. The junctions were fabricated by electron-beam lithography using 
a lift-off technique. Pt films, 20 nm thick, were deposited on a Si/SiC>2 substrate by 
magnetron sputtering with approximately a 15° deposition angle with respect to the 
surface normal of the sample. The samples were rotated about the surface normal 
during deposition. The Pt deposition rate was 0.028 nm/s and the base pressure was 
5 x 10~ 8 Torr; an Ar pressure of 1.5 mTorr was used during deposition. The junction 
thickness at its narrowest point was approximately 8 nm to 10 nm because of a shadow 
effect during the sputtering deposition process, as shown in figure [H(a). The junction 
with this kind bow-tie shape tends to break only in the middle, does not require high 
voltages to break [figures [T^b) and (c)], and can be broken by an electrical pulse with a 
low negative voltage, as shown in figure [H(d) [18J. The total resistance of the junctions 
(including parasitic resistance) at room temperature ranged between 100 Q and 2000 Q. 
The narrowest gap between electrodes formed in this way are on the order of 5 nm, as 
shown in figure [Q 
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Figure 1. (a) Atomic force microscopy (AFM) image for one Pt junction which 20 
nm thickness in large area, 8 to 10 nm in the middle, (b) AFM image after breaking 
the junction by electromigration at room temperature, (c) I — V measurement for 
breaking the junction at 0.9 V and room temperature, (d) AFM image for another 
junction broken by an electrical pulse (-50 mV) at T = 6 K. 



2.2. Protein incorporation and viability 

Mb and apoMb at least 95% pure were obtained from horse skeletal muscle (Sigma) 
in its powder state and essentially salt-free. The protein powders were dissolved in a 
10 niM phosphate buffer saline (PBS, Sigma, pH = 7.4) with 14.2 mM and 31.5 fiM 
concentrations for Mb and apoMb, respectively. Protein solutions were passed through 
a 0.2 /im pore size filter to remove larger protein clusters that may have formed. In 
order to remove possible contaminants from the electrode surface, junctions were rinsed 
with optima acetone and isopropanol and subsequently cleaned with O2 plasma. After 
the cleaning procedure, the protein solution was dropped on the samples. The samples 
were kept at 4 °C in order to let water evaporate. A drop of the solution was placed on 
the unbroken junction prior to cooling down the samples. 

Linker molecules were avoided because junctions exposed to them (without protein) 
generated resonant tunneling electronic data. These data showed a weak dependence on 
gate voltage, which would have made it difficult to disntiguish the protein behavior from 
that of the linker (see supplementary information). Raman spectra showed that the Mb 
was active after cooling down to T = 5 K and then warming up to room temperature [19J. 
Therefore, it is unlikely that the protein denatured during the cooling process. Another 
possible problem is that the local temperature during the electromigration process might 
be high enough to denature the protein [20, ED, 122] ■ Measurements of the electrical 
resistance at low voltages as a function of temperature yielded an estimate of a maximum 
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Figure 2. Sketch of three-terminal device and AFM image of a bare Pt junction 
broken by electromigration. The AFM image indicates that the gap is on the order of 
5 nm wide. 

temperature of approximately 400 K, which is consistent with estimates of 495 K for 
junctions broken at room temperature [23]. Additional Raman scattering measurements 
performed on protein heated up in air in a furnace or hot plate demonstrated that 
myoglobin is viable after heating to 400 K (see supplementary information). 

The junctions with a bow-tie configuration were fabricated from Pt thin films, 
8 to 10 nm thick in the middle, 200 nm wide and deposited on a Si wafer with 
a 200 nm thick Si0 2 insulating layer. The doped Si substrate was used as a gate 
contact in all samples, as shown in figure [2j Electrical measurements were performed 
immediately after breaking the junction at low temperature (either 77 K or 6 K). 
The conductance (dl/dV) as a function of V# and Vq was acquired by numerically 
differentiating current-bias voltage {I-V) curves. Junctions were measured without any 
protein, with apomyoglobin, and with myoglobin. 

3. Results 

3.1. Bare junctions 

In the control experiments with bare junctions (without protein), 107 Pt junctions 
were broken by electromigration at T = 77 K. The broken junctions had a resistance 
ranging from 1 MQ to 1 GQ. All but one of the broken junctions did not show 
differential conductance peaks. Figure |3]^a) is a typical result measured at T = 77 K. 
The conductivity was relatively low and there was a significant amount of noise for 
these structures. Only one sample had conductance peaks, as shown in figure Mjo)- 
In addition, 31 bare Pt junctions were broken by electromigration and measured at 
T = 5.5 K. Seventeen of these did not have conductance peaks as shown in figure M^c). 
Out of the remaining fourteen junctions that had conductance peaks, eleven broken 
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Figure 3. Typical results for bare Pt junctions, (a) Conductance as a function of 
bias voltage measured at T = 77 K. No conductance peaks are found near zero bias, 
(b) Conductance of the only sample out of 107 junctions broken at T — 77 K that had 
noticeable conductance features, (c) Typical differential conductance at T = 5.5 K for 
bare junctions. No conductance peaks are evident, (d) SET behavior measured in a 
bare Pt junction at T = 5 K [gray scale from -3.5 x 1CT 7 S (black) to 1.8 x 1CT 5 S 
(white)]. 

junctions did not have a dependence on the gate voltage. The remaining three junctions 
did have conductance peaks whose positions were dependent on the gate voltage. 
Coulomb diamonds, characteristic of SET behavior, were observed in the conductance 
as a function of gate and bias voltages as shown in figure |3]^d). This demonstrates that 
occasionally Pt grains are formed, as has been previously observed in Au junctions, 
which can mask the molecular information |16[ [HI [22] . These islands have several 
energy levels and the conductance peaks from metal grains can be tuned by gate 
voltage modulation. Many Coulomb diamonds were observed, each corresponding to 
the addition of an electron to the island. It is important to note that the appearance 
of several Coulomb diamonds is unlikely to happen in single-protein junctions, since 
excessive charging can denature the protein. 

3.2. Apomyoglobin experiments 

In the apoMb experiments, out of 125 Pt protein-coated junctions broken at T = 77 K, 
nine had conductance peaks with a very weak response to Vq. Figure H^a) is the 
conductance measured at T = 77 K for one of these apoMb samples. One conductance 
peak was approximately constant as a function of Vq. This peak, measured at Vq = as 
a function of temperature, is shown in figure Hb- Figure H^c) shows how the peak height 
depended on temperature. The monotonic decrease of the peak height with increasing 
temperature, together with saturation at low temperatures (T < 100 K) indicates that 
tunneling, and not Mott hopping, is the dominant electron conduction mechanism in 
this device [211 [25]. Tunneling is observed because the conductance is smaller than the 
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conductance quantum of G a = 2e 2 /h = 7.75 x 10~ 5 S. The lack of dependence of Vq 
means, however, that the separation of energy levels inside the protein are much larger 
than k B T [26]. 

Forty-three apoMb samples were measured at T R3 6 K. Fifteen of them had 
conductance peaks and two had a very weak V G dependence. For three of the apoMb 
samples, conductance lines changed slightly with V G as shown in figure IU(d). The 
spectra measured at V G = at several temperatures is shown in figure 11(e). The zero 
conductance region near Vb = may arise from a Coulomb blockade effect, but the lack 
of Coulomb triangles indicates that the energy levels in the apoMb protein are too close 
in energy to be able to observe a Coulomb staircase. The temperature dependence of the 
the conductance at Vb — is shown in figure H](f). For temperatures below T = 20 K the 
conductance is independent of temperature, characteristic of a tunneling process, but at 
higher temperatures the conductance quickly increases until the spectrum completely 
loses its tunneling gap above T = 100 K. The line that starts at V B ~ 132.5 mV, 
indicated by the red arrow, may arise from the vibrational spectra of amino acid redox 
centers [2"7] . 

3.3. Myoglobin experiments 

A total of 215 Mb samples were measured at T rj 5 K, out of which 21 had conductance 
peaks that strongly depended on V G . In another approach, junctions were broken by 
a single electrical pulse with a very low negative voltage with a mechanism similar to 
voltage-pulse-induced electromigration [18]. These thin (~ 8 nm) Pt junctions can be 
broken by a low negative voltage « —60 mV which avoids potentially heating the sample 
to high temperatures. Out of 64 Pt Mb-coated junctions broken by electrical pulse at 
T = 6 K, four had an obvious dependence on Vq. 

Of the 25 samples that had gate voltage-dependent behavior, the conductance maps 
can be divided into three types: 1) data where a Kondo-like resonance is observed, that 
is, conductance near zero bias voltage that is independent of gate voltage; 2) Coulomb 
blockade without a two-step process; and 3) Coulomb blockade with a two-step process. 

Figure [5] shows results from one of the samples with a Kondo-like resonance. When 
the junction was initially broken by a negative electrical pulse, no conductance peaks 
were found. However, there was a sudden increase in the current when the voltage 
was increased again, indicating a modification of the contact and possible protein 
incorporation in the gap of electrodes [ID] , although this could also be due to Joule 
heating resulting from a sudden narrowing of the junction in the bowtie region [28J. 
Conductance peaks were observed thereafter. In the differential conductance map, 
shown in figure [5](a), two intense peaks around Vb — are independent of V G for 
< Vg < 3 V. For higher Vq values, the peaks are composed of half a Coulomb blockade- 
like triangle. This could be an indication of a Kondo-resonance assisted tunneling [10J, 
although verification would require measurement in large magnetic fields. The two 
lines indicated by the green arrows, corresponding to green arrows in figure [5](a), do 
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Figure 4. ApoMb results. (a) Differential conductance (dl/dV) measured at 
T = 77 K as a function of bias and gate voltages. Gray scale from 1.56 x 10 -5 S (black) 
to 2.64 x 10 -5 S (white) . (b) dl/dV measured at Vq = for various temperatures, 
(c) dl/dV of peak in (b) as a function of 1/T. (d) dl/dV for another Mb sample [gray 
scale from (black) to 5.0 x 10~ 6 S (white)] measured at T = 6 K. The red arrow 
points to a possible vibration-assisted conduction line, (e) dl/dV at Vg = for various 
temperatures, (f) Conduction at Vb = and Vg = for sample in (d) as a function 
of 1/T. 



not depend strongly on V G for < Vq < 1 V. The conductance feature starting at 
Vb ~ 30 mV could arise from vibrational excitations as this energy coincides with 
the energy of the Fe-His vibration (27.3 meV). Conductance steps are also present as 
indicated in figure[5](b). These steps result in peaks in d 2 I/dV 2 , indicated by the black (9 
mV), blue (35 mV), and magenta (~ 80 mV) arrows in figure 0(b). Taking into account 
that modes increase in energy at low T, the arrows correspond to heme- doming [29, 30J, 
Fe-His, and, heme-group deformation (83.7 meV) modes, respectively [31] . 

Conductance data for a sample which displayed Coulomb blockade behavior without 
a two-step process are shown in figure [6](a). Figure Mjo) is a calculation using reasonable 
parameters that reproduce the main characteristics of the experimental data. The 
appearance of two wide conductance lines starting at Vb ~ 5.2 mV, indicated by red 
arrows in figure [6](a), arise from inelastic tunneling with the 4.96 meV vibrational mode 
of the heme group [291 E0] ■ 

Another sample with SET behavior is shown in figure [6](b). In this case, two 
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Figure 5. Mb results, (a) dl/dV as function of Vb and Vg at T = 6 K [gray scale is 
(black) to 1.33 x 10 -5 S (white)]. Green arrows point to conduction lines corresponding 
to Fe-His vibration-assisted tunneling, (b) dl/dV (black) and d 2 I/dV 2 (red) spectra 
at Vg = 11 V. Vertical arrows correspond to inelastic tunneling peaks. 

deformed Coulomb blockade triangles near zero bias voltage can be found without a 
degeneracy point in the differential conductance graph. This means that two charge 
states do not have the same energy, and in this case the separation between the top 
triangles is > 60 mV. This separation, as shown by the calculation in figure Etc), is an 
indication that the two-step electron transfer process is a viable mechanism. In addition, 
there are three areas that exhibit negative differential conductance (NDC), indicated 
by the three blue arrows. One explanation of how NDC behavior may arise is from 
the existence of two quantum dots in series . This suggests that these measurements 
were the result of transport through two or more myoglobin proteins in series with each 
other. The faint line at Vb ~ —25 mV, indicated by the green arrow, does not have a 
strong dependence on Vg, which indicates electron tunneling by the medium around the 
heme group. We also note that the conductance feature that starts at Vb = —97 mV, 
indicated by the red arrow, could arise from one of the excited electronic level of the 
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Figure 6. (a) Differential conductance (dl/dV) data from a Mb sample [gray scale 
from (black) to 5.0 x 1CP 6 S (white)] at T = 6 K. The red arrows point to vibration- 
assisted conduction lines, (b) Differential conductance T = 6 K [gray scale from 
-0.27 x 10~ 6 S (black) to 1.6 x 10~ 6 S (white)] for another Mb sample. The blue 
arrows point to regions of negative differential conductance (NDC). The green arrow 
points to a faint conductance line independent of Vg ■ The red arrow points to another 
conductance line possibly arising from a vibrational-assisted tunneling process or an 
excited electronic level in the protein(s). 



protein or a vibrational excitation of the protein. In particular, this energy coincides 
with the energy of the tryptophan amino acid in Mb which is known to have an energy 
of 760 cm -1 = 94.2 mV (measured at room temperature) [32]. The wide conductance 
line that starts at Vb ~ 93 mV, indicated by the pink arrow, is another line that could 
arise from the excited state or vibrational excitation. 
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Figure 7. Heme group structure (a) and energies of Fe d-orbitals in dry metMb (b). 
4. Discussion 

Here we discuss our results in terms of the known properties of myoglobin. Figure [7(a) 
shows the heme group structure of dry Mb. The fifth Fe ligand binds covalently with 
amino acid His 93, and therefore the Fe-His vibrational mode (27.3 meV) is an excellent 
heme-group indicator [29]. The sixth Fe ligand is occupied by H 2 for aqueous Mb but 
is coordinated with His 64 when dry. In the latter case the Fe atom is in the plane of 
porphyrin ring in a low spin (S = 1/2) configuration (metMb) [31]. The ligand field 
splits the Fe energy levels into two e g and three t2 g orbitals, as illustrated in figure [7(b), 
with an energy gap A « 1.0 eV for metMb [29]. These orbitals split further by 0.25 eV 
to 0.37 meV due to molecular orbital interactions [31], which is difficult to detect even at 
a cryogenic temperature of T ~ 5 K. Excited states are also difficult to probe because of 
low Vb (< 300 mV) used in the measurements (higher voltages destroyed our samples). 

A Mb SET can be modeled as a single energy level device with inelastic tunneling 
processes possible due to the localized vibrational modes of the heme group. However, 
it is known that redox reactions can cause significant conformational changes of the 
heme group [33] and that the associated relaxation time could range from picoseconds 
to seconds [34J. If the relaxation time is long, the electron resonantly tunnels while 
the protein remains in the intermediate state and its degeneracy point is at Vb — 0. 
If the relaxation time is short (on the order of a few ps), the protein quickly relaxes 
from the intermediate state to its ground state, whose level may fall below the Fermi 
level of electrodes, as illustrated in figure Eta). In this case, the electron is trapped 
in the protein, current flow is blocked, and the Coulomb blockade lacks a degeneracy 
point. This is further illustrated by the calculations in figures [E(b) and [8(c) performed by 
solving tunneling rate equations [351 [36] (see Appendix A). When the two-step relaxation 
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Figure 8. (a) Two-step electron transfer by redox center. SE is the difference in 
conformational energies between the unoccupied and occupied states. The Fermi levels 
of the source and drain electrodes are indicated by /is and /id, respectively, (b) 
Calculated differential conductance dl/dV as a function of bias and gate voltages Vb 
and V G for T = 6 K, SE = mV, V c = 0.307 V, C D : C s : C G = 24 : 40 : 1, tunneling 
rate 10 GHz (between the protein and the source or drain electrodes). Cd, Cs and 
Cg are the capacitances of the drain, source, and gate electrical contacts, respectively. 
There is no two-step process in this case (SE = 0). (c) Calculated dl/dV with a two- 
step process for SE = 10 mV, other parameters same as in (b). See Appendix A for 
description of calculation. 



process does not take place, normal Coulomb triangles are observed in the differential 
conductance dl/dV map. When the relaxation process occurs, a gap opens at the 
intersection of the Coulomb triangles (degeneracy point). 
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It is important to note that amino acids around the heme group may complicate the 
analysis by acting as intermediaries for electron hopping processes and redox centers [27J . 
These aromatic amino acids have low frequency vibrational modes, such as phenylalanine 
(124 meV), tyrosine (79.9 meV, 103 meV, and 106 meV), and tryptophan (94.2 meV, 
109 meV and 125 meV) [32]. 

A fundamental question in analyzing our results is determining whether the protein 
is indeed being measured. Regarding the results of the bare junctions with respect to 
the Mb junctions, we notice that the conductivity maps in figures [5] and [6] do not show 
multiple degeneracy point like those of the bare junctions in figure E^d), and the energy 
scale (gate and bias voltages) used are completely different. This is an indication that 
the myoglobin data come from the protein and not from metallic islands formed in 
the junctions. We note that a similar observation has been made for break junctions 
used to measure single OPV-5 molecules, where gold grains remaining in the junction 
create multiple Coulomb diamonds which are superimposed on the signal from the 
molecule p2]. As to whether the heme group plays a role in resonant tunneling processes, 
we note that none of the apoMb samples' conductance showed a strong dependence on 
gate voltage, while only Mb samples had this dependence. Another indication that the 
data represent the conductance of the protein are the inelastic tunneling features of the 
conductance plots in figures [5] and El As mentioned above, these features correspond to 
the known vibrational energies of the heme group. 

The effects of the protein orientation can be surmised from inelastic vibrational 
spectra. Two pathways for electron transfer between the heme iron and electrodes are 
through the porphyrin ring or the His93 group |31j. Transfer through the porphyrin 
ring occurs as a result of porphyrin rings p-orbital hybridization with the ci-orbital of 
iron (31] combined with a slow porphyrin- iron structural relaxation (37]. In this case, the 
resonant energy level does not change during tunneling, and results like those in figures 
(8](b) and El(a) are expected. On the other hand, the dip near Vb = and Vq < 3 V in 
Fig. Mb) could be due to the very fast relaxation ( ps) of the Fe-His bond (37J [38] - The 
qualitative difference between figures Eta) and Eta) could be due to different orientations 
of the heme group with respect to the electrodes. Moreover, because the heme group 
is not at the center of the protein [15], certain orientations may cause the heme to be 
closer to the electrodes, thus increasing the probability of observing Kondo resonances. 

Another complication is that the large lateral size of the break junctions could result 
in several proteins with different conformations being measured in parallel, or in protein 
dimers consisting of two proteins measured in series if the gap is somewhat larger than 
a single protein. We assume that conductance is dominated by the narrowest portions 
of the break junction since tunneling processes usually depend exponentially on the 
effective width of the barriers, and therefore it is unlikely that trimers or larger protein 
agglomerations will play much of a role in the conductance data. The data shown in 
figure E^a) agree well with what is expected from a single protein conductance, whereas 
the data in figure El^b) show indications of multiple protein conductance (either in series 
or in parallel), in addition to a two-step tunneling process. Analyzing the data with 
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more sophisticated models that take into account multiple protein conductance may be 
helpful in further understanding the conductance mechanism, but this is beyond the 
scope of this paper. 

5. Conclusions 

We have presented results from a series of experiments on metallic break junctions 
with apomyoglobin and myoglobin. Some myoglobin junctions displayed single-electron 
transistor behavior consistent with resonant tunneling through the heme group. In one 
case, tunneling appears to be a result of a two-step process where the protein changes 
conformation during the tunneling process. Elastic and intelastic tunneling data are also 
consistent with resonant tunneling from known vibrational modes of the protein. Our 
results underscore the importance of tunneling, as opposed to hopping, as an electron 
transfer mechanism in redox proteins. 

Future challenges include determining the spin state of a single metalloprotein via 
a SET configuration for sensor applications and using other metalloproteins, such as 
cytochrome P450, that have binding properties to specific molecules. This will require a 
greater yield of successful devices and control of the protein orientation and conformation 
on the electrodes. 
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Appendix A. Rate Equations for Two-step Electron Transfer 

In the following treatment we use the same approach as in references [35] and [36]. The 
rate equation for the transition of an electron into an energy level in a SET is 

dP a 

— ^(Tp^aPp — Fa^pPa), (A.l) 

where P is the probability that the molecule is in a specific state and and are 
the transition rates between a and (3 states. The current can be calculated by electron 
transfer through either of the two tunnel barriers corresponding to the source and drain 
electrical contacts. The two currents at equilibrium are equal. The current through the 
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source is 



(A.2) 



where e is the charge of the electron, T^gP a is the contribution of the source to T a _+p 
multiplied by ±1 depending on the transition contribution to the current. 

Figure [8(a) is a diagram representing the two-step resonant tunneling process. The 
proteins ferric state (Fe +3 ) is labeled by (0) and the ferrous state (Fe +2 ) is labeled 
(1). Ts and are the tunneling rates between the protein and the source and drain 
electrodes, respectively. The energy difference between the ferric and ferrous states is 
SE. We define Pq and Pi as the probabilities that the level in the quantum dot (protein) 
are unoccupied and occupied, respectively. 

For single electron transport, the total tunneling time r can be estimated from 
t = e/I, so that r ~ps for a current of I ~ 10~ 7 A. We therefore set up the rate 
equations for two-step electron transfer by making the following assumptions: 

(i) When the electron tunnels into the protein, the energy level relaxes to level 1 via 
a structural distortion in a time ~ r. 

(ii) The process is not reversible and fluctuations between levels 1 and resulting 
from nuclear distortions are neglected. This is true at low temperatures where 
the thermal energy is much smaller than the vibrational energy of the heme group 
(fc B T«£ vib ). 

With these assumptions, the rate of change of Pq and P\ can be written as 



dPo 
dt 



-p (v s f° s + r D f D )+p 1 [t s (i - f s ) +r D (i- f D 



dPi 



(A.3) 



and f D are the Fermi distribution functions of the source and drain 



where fg 

electrodes, respectively, for the unoccupied (occupied) energy levels. These can be 
written as 

-l 



Js. 



S,D 



fs,D 



1 + exp 



1 + exp 



(/' 



N- 



_i - 5E - /i SjD ) 



k B T 



(A.4) 



(A.5) 



Here fi^ is the chemical potential of the system with N electrons and /ijv+i is the 
chemical potential when an additional electron is added to the system. The latter can 
be expressed as 

CqVg + CsVsd 



fJ'N+i = Eq — eVj ot = E — e- 



a 



(A.6) 



total 



where V^ Q ^ is the potential of the quantum dot, Vsd is the source-drain potential 
difference, Vq is the gate potential, C to tai = Co + Cs + Cd is the total capacitance 
of th all the electrodes in the system with Co, Cs, and Cd being the capacitances of 
the gate, source, and drain electrodes. The energy E = eVcCc/Ctotai is proportional 
to the crossing potential Vc, which is the potential at which current can flow at low bias 
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due to matching of the empty electron energy level with the source electrode chemical 
potential. 

In equilibrium, = = and Pq + P\ = 1. Combining this equilibrium 
condition with equation IA.9I and the expression for the current 



leads to 



and 



/ = e (P 1 F S (1 - fs) - PoTsfs) (A.7) 



r s (i - /j) + r c (i - /£) 
r s (i + /°-/i) + r D (i + /o-/i) 1 f 



er, 



(A - 9) 

In figures EJ^b) and (c) dl/dV was calculated by numerically taking the derivative 
of the expression for J in equation I A. 9 1 with respect to V$d- As in the experiments, the 
chemical potential of the drain was set to fio = (i.e., the drain contact was grounded). 
Note that if 5E = there is no two-step process and the result of a single-level quantum 
dot is recovered. 
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